Differential amplification of the subtelomeric satellite DNA JcSAT1 in the genus Jatropha L. (Euphorbiaceae).

Satellite DNAs (satDNAs) are highly repetitive sequences that occur in virtually all eukaryotic genomes and can undergo rapid copy number and nucleotide sequence variation among relatives. After chromosomal mapping of the satDNA JcSAT1, it was found a large accumulation at subtelomeres of Jatropha curcas (subgenus Curcas), but an absence of these monomers in J. integerrima (subgenus Jatropha). This fact suggests a dynamic scenario for this satellite repeat in Jatropha genomes. Here, we used a multitasking approach (sequence analysis, DNA blotting and chromosomal mapping) to investigate the molecular organization and chromosomal abundance and distribution of JcSAT1 in a broader group of species from the subgenus Jatropha (J. gossypiifolia, J. mollissima, J. podagrica, and J. multifida) in addition to J. curcas, with the aiming of understanding the evolution of this satDNA. Based on the analysis of BAC clone sequences of J. curcas, a large array (~ 30 kb) of 80 homogeneous monomers of JcSAT1 was identified in BAC 23J11. The monomer size was conserved (~ 358 bp) and contained a telomeric motif at the 5' end. PCR amplification coupled with a Southern blot revealed the presence of JcSAT1-like sequences in all species examined. However, a large set of genome copies was identified only in J. curcas, where a ladder-like pattern with multimers of different sizes was observed. In situ hybridization of BAC 23J11 confirmed the subtelomeric pattern for J. curcas, but showed no signals on chromosomes of species from the subgenus Jatropha. Our data indicate that JcSAT1 is a highly homogeneous satDNA that originated from a region near the telomeres and spread throughout the chromosomal subtermini, possibly due to frequent ectopic recombination between these regions. The abundance of JcSAT1 in the genome of J. curcas suggests that an amplification event occurred either at the base of the subgenus Curcas or at least in this species, although the repeat is shared by all species of the genus studied so far.

Evolvability and constraint in the evolution of three-dimensional flower morphology.

PREMISE: Flower phenotypes evolve to attract pollinators and to ensure efficient pollen transfer to and from the bodies of pollinators or, in self-compatible bisexual flowers, between anthers and stigmas. If functionally interacting traits are genetically correlated, response to selection may be subject to genetic constraints. Genetic constraints can be assessed by quantifying standing genetic variation in (multivariate) phenotypic traits and by asking how much the available variation is reduced under specific assumptions about phenotypic selection on functionally interacting and genetically correlated traits. METHODS: We evaluated multivariate evolvability and potential genetic constraints underlying the evolution of the three-dimensional structure of Dalechampia blossoms. First, we used data from a greenhouse crossing design to estimate the G matrix for traits representing the relative positions of male and female sexual organs (anthers and stigmas) and used the G matrix to ask how genetic variation is distributed in multivariate space. To assess the evolutionary importance of genetic constraints, we related standing genetic variation across phenotypic space to evolutionary divergence of population and species in the same phenotypic directions. RESULTS: Evolvabilities varied substantially across phenotype space, suggesting that certain traits or trait combinations may be subject to strong genetic constraint. Traits involved functionally in flower-pollinator fit and autonomous selfing exhibited considerable independent evolutionary potential, but population and species divergence tended to occur in phenotypic directions associated with greater-than-average evolvability. CONCLUSIONS: These results are consistent with the hypothesis that genetic constraints can hamper joint trait evolution towards optimum flower-pollinator fit and optimum self-pollination rates.

Evolution and function of ubiquitin-specific proteases (UBPs): Insight into seed development roles in plants.

The tung oil produced by the tung tree (Vernicia fordii) provides resources for the manufacture of biodiesel. Ubiquitin-specific proteases (UBPs) are the largest group of deubiquitinases and play key roles in regulating development and stress responses. Here, 21 UBPs were identified in V. fordii, roughly one-half the number found in Manihot esculenta and Hevea brasiliensis. Most UBP duplications are produced from whole-genome duplication (WGD), and significant differences in gene retention existed among Euphorbiaceae. The great majority of UBP-containing blocks in V. fordii, V. montana, Ricinus communis, and Jatropha curcas exhibited extensive conservation with the duplicated regions of M. esculenta and H. brasiliensis. These blocks formed 14 orthologous groups, indicating they shared WGD with UBPs in M. esculenta and H. brasiliensis, but most of these UBPs copies were lost. The UBP orthologs contained significant functional divergence which explained the susceptibility of V. fordii to Fusarium wilt and the resistance of V. montana to Fusarium wilt. The expression patterns and experiments suggested that Vf03G1417 could affect the seed-related traits and positively regulate the seed oil accumulation. This study provided important insights into the evolution of UBPs in Euphorbiaceae and identified important candidate VfUBPs for marker-assisted breeding in V. fordii.

Identification of Gene Biomarkers for Tigilanol Tiglate Content in Fontainea picrosperma.

Tigilanol tiglate (EBC-46) is a small-molecule natural product under development for the treatment of cancers in humans and companion animals. The drug is currently produced by purification from the Australian rainforest tree Fontainea picrosperma (Euphorbiaceae). As part of a selective-breeding program to increase EBC-46 yield from F. picrosperma plantations, we investigated potential gene biomarkers associated with biosynthesis of EBC-46. Initially, we identified individual plants that were either high (>0.039%) or low EBC-46 (<0.008%) producers, then assessed their differentially expressed genes within the leaves and roots of these two groups by quantitative RNA sequencing. Compared to low EBC-46 producers, high-EBC-46-producing plants were found to have 145 upregulated genes and 101 downregulated genes in leaves and 53 upregulated genes and 82 downregulated genes in roots. Most of these genes were functionally associated with defence, transport, and biosynthesis. Genes identified as expressed exclusively in either the high or low EBC-46-producing plants were further validated by quantitative PCR, showing that cytochrome P450 94C1 in leaves and early response dehydration 7.1 and 2-alkenal reductase in roots were consistently and significantly upregulated in high-EBC-46 producers. In summary, this study has identified biomarker genes that may be used in the selective breeding of F. picrosperma.

Dioecy and chromosomal sex determination are maintained through allopolyploid speciation in the plant genus Mercurialis.

Polyploidization may precipitate dramatic changes to the genome, including chromosome rearrangements, gene loss, and changes in gene expression. In dioecious plants, the sex-determining mechanism may also be disrupted by polyploidization, with the potential evolution of hermaphroditism. However, while dioecy appears to have persisted through a ploidy transition in some species, it is unknown whether the newly formed polyploid maintained its sex-determining system uninterrupted, or whether dioecy re-evolved after a period of hermaphroditism. Here, we develop a bioinformatic pipeline using RNA-sequencing data from natural populations to demonstrate that the allopolyploid plant Mercurialis canariensis directly inherited its sex-determining region from one of its diploid progenitor species, M. annua, and likely remained dioecious through the transition. The sex-determining region of M. canariensis is smaller than that of its diploid progenitor, suggesting that the non-recombining region of M. annua expanded subsequent to the polyploid origin of M. canariensis. Homeologous pairs show partial sexual subfunctionalization. We discuss the possibility that gene duplicates created by polyploidization might contribute to resolving sexual antagonism.

Overexpression of PvFAD3 Gene from Plukenetia volubilis Promotes the Biosynthesis of α-Linolenic Acid in Transgenic Tobacco Seeds.

The ω-3 fatty acid desaturase (FAD3) gene encodes a rate-limiting enzyme in the synthesis of α-linolenic acid. In this study, homologous cloning was used to obtain the full-length sequence of the PvFAD3 gene of Plukenetia volubilis. The full-length DNA sequence was 1871 bp long, with 8 exons and 7 introns. The structural analysis of the amino acid sequence revealed that the PvFAD3 protein contained three histidine-conserved regions and an endoplasmic reticulum retention signal. The real-time reverse transcription-polymerase chain reaction performed for determining the expression patterns of the PvFAD3 gene in different tissues of P. volubilis showed that PvFAD3 expression was highly expressed in the fast oil accumulation stage of seed. The analysis of subcellular localization assay in epidermal cells of tobacco (Nicotiana benthamiana) leaves showed that the PvFAD3 protein was mainly localized in the endoplasmic reticulum. Seed-specific overexpression vectors were constructed, and Agrobacterium-mediated genetic transformation was performed to obtain transgenic tobacco plants overexpressing PvFAD3. The results of fatty acid assays performed using harvested seeds showed a significant increase in α-linolenic acid content, a dramatic decrease in linoleic acid content, and an obvious increase in oil content in transgenic tobacco seeds. Collectively, the PvFAD3 gene of P. volubilis was confirmed as a key enzyme gene for α-linolenic acid synthesis; thus, indicating that the PvFAD3 gene can be used for fatty acid fraction improvement in oilseed plants.

Recurrent allopolyploidization, Y-chromosome introgression and the evolution of sexual systems in the plant genus Mercurialis.

The plant genus Mercurialis includes dioecious, monoecious and androdioecious species (where males coexist with hermaphrodites). Its diversification involved reticulate evolution via hybridization and polyploidization. The Y chromosome of the diploid species Mercurialis annua shows only mild signs of degeneration. We used sequence variation at a Y-linked locus in several species and at multiple autosomal and pseudoautosomal loci to investigate the origin and evolution of the Y chromosome across the genus. Our study provides evidence for further cases of allopolyploid speciation. It also reveals that all lineages with separate sexes (with one possible exception) share the same ancestral Y chromosome. Surprisingly, males in androdioecious populations of hexaploid M. annua carry a Y chromosome that is not derived from either of its two putative progenitor lineages but from a more distantly related perennial dioecious lineage via introgression. These results throw new light on the evolution of sexual systems and polyploidy in Mercurialis and secure it as a promising model for further study of plant sex chromosomes. This article is part of the theme issue 'Sex determination and sex chromosome evolution in land plants'.

Genomic insights into lineage-specific evolution of the oleosin family in Euphorbiaceae.

BACKGROUND: Lipid droplets (LDs) present in land plants serve as an essential energy and carbon reserve for seed germination and seedling development. Oleosins, the most abundant structural proteins of LDs, comprise a small family involved in LD formation, stabilization and degradation. Despite their importance, our knowledge on oleosins is still poor in Euphorbiaceae, a large plant family that contains several important oil-bearing species. RESULTS: To uncover lineage-specific evolution of oleosin genes in Euphorbiaceae, in this study, we performed a genome-wide identification and comprehensive comparison of the oleosin family in Euphorbiaceae species with available genome sequences, i.e. castor bean (Ricinus communis), physic nut (Jatropha curcas), tung tree (Vernicia fordii), Mercurialis annua, cassava (Manihot esculenta) and rubber tree (Hevea brasiliensis), and a number of five, five, five, five, eight and eight members were found, respectively. Synteny analysis revealed one-to-one collinear relationship of oleosin genes between the former four (i.e. castor bean, physic nut, tung tree and M. annua) as well as latter two species (i.e. cassava and rubber tree), whereas one-to-one and one-to-two collinear relationships were observed between physic nut and cassava, reflecting the occurrence of one recent whole-genome duplication (WGD) in the last common ancestor of cassava and rubber tree. The presence of five ortholog groups representing three previously defined clades (i.e. U, SL and SH) dates back at least to the Malpighiales ancestor, because they are also conserved in poplar (Populus trichocarpa), a tree having experienced one Salicaceae-specific recent WGD. As observed in poplar, WGD was shown to be the main driver for the family expansion in both cassava and rubber tree. Nevertheless, same retention patterns of WGD-derived duplicates observed in cassava and rubber tree are somewhat different from that of poplar, though certain homologous fragments are still present in rubber tree. Further transcriptional profiling revealed an apparent seed-predominant expression pattern of oleosin genes in physic nut, castor bean and rubber tree. Moreover, structure and expression divergence of paralogous pairs were also observed in both cassava and rubber tree. CONCLUSION: Comparative genomics analysis of oleosin genes reported in this study improved our knowledge on lineage-specific family evolution in Euphorbiaceae, which also provides valuable information for further functional analysis and utilization of key members and their promoters.

Triacylglycerol biosynthesis in shaded seeds of tung tree (Vernicia fordii) is regulated in part by Homeodomain Leucine Zipper 21.

Light quantity and quality affect many aspects of plant growth and development. However, few reports have addressed the molecular connections between seed oil accumulation and light conditions, especially dense shade. Shade-avoiding plants can redirect plant resources into extension growth at the expense of leaf and root expansion in an attempt to reach areas containing richer light. Here, we report that tung tree seed oil accumulation is suppressed by dense shade during the rapid oil accumulation phase. Transcriptome analysis confirmed that oil accumulation suppression due to dense shade was attributed to reduced expression of fatty acid and triacylglycerol biosynthesis-related genes. Through weighted gene co-expression network analysis, we identified 32 core transcription factors (TFs) specifically upregulated in densely shaded seeds during the rapid oil accumulation period. Among these, VfHB21, a class I homeodomain leucine zipper TF, was shown to suppress expression of FAD2 and FADX, two key genes related to α-eleostearic acid, by directly binding to HD-ZIP I/II motifs in their respective promoter regions. VfHB21 also binds to similar motifs in the promoters of VfWRI1 and VfDGAT2, two additional key seed lipid regulatory/biosynthetic genes. Functional conservation of HB21 during plant evolution was demonstrated by the fact that AtWRI1, AtSAD1, and AtFAD2 were downregulated in VfHB21-overexpressor lines of transgenic Arabidopsis, with concomitant seed oil reduction, and the fact that AtHB21 expression also was induced by shade. This study reveals some of the regulatory mechanisms that specifically control tung tree seed oil biosynthesis and more broadly regulate plant storage carbon partitioning in response to dense shade conditions.

The P450 multigene family of Fontainea and insights into diterpenoid synthesis.

BACKGROUND: Cytochrome P450s (P450s) are enzymes that play critical roles in the biosynthesis of physiologically important compounds across all organisms. Although they have been characterised in a large number of plant species, no information relating to these enzymes are available from the genus Fontainea (family Euphorbiaceae). Fontainea is significant as the genus includes species that produce medicinally significant epoxy-tigliane natural products, one of which has been approved as an anti-cancer therapeutic. RESULTS: A comparative species leaf metabolome analysis showed that Fontainea species possess a chemical profile different from various other plant species. The diversity and expression profiles of Fontainea P450s were investigated from leaf and root tissue. A total of 103 and 123 full-length P450 genes in Fontainea picrosperma and Fontainea venosa, respectively (and a further 127/125 partial-length) that were phylogenetically classified into clans, families and subfamilies. The majority of P450 identified are most active within root tissue (66.2% F. picrosperma, 65.0% F. venosa). Representatives within the CYP71D and CYP726A were identified in Fontainea that are excellent candidates for diterpenoid synthesis, of which CYP726A1, CYP726A2 and CYP71D1 appear to be exclusive to Fontainea species and were significantly more highly expressed in root tissue compared to leaf tissue. CONCLUSION: This study presents a comprehensive overview of the P450 gene family in Fontainea that may provide important insights into the biosynthesis of the medicinally significant epoxy-tigliane diterpenes found within the genus.

Using ecological context to interpret spatiotemporal variation in natural selection.

Spatiotemporal variation in natural selection is expected, but difficult to estimate. Pollinator-mediated selection on floral traits provides a good system for understanding and linking variation in selection to differences in ecological context. We studied pollinator-mediated selection in five populations of Dalechampia scandens (Euphorbiaceae) in Costa Rica and Mexico. Using a nonlinear path-analytical approach, we assessed several functional components of selection, and linked variation in pollinator-mediated selection across time and space to variation in pollinator assemblages. After correcting for estimation error, we detected moderate variation in net selection on two out of four blossom traits. Both the opportunity for selection and the mean strength of selection decreased with increasing reliability of cross-pollination. Selection for pollinator attraction was consistently positive and stronger on advertisement than reward traits. Selection on traits affecting pollen transfer from the pollinator to the stigmas was strong only when cross-pollination was unreliable and there was a mismatch between pollinator and blossom size. These results illustrate how consideration of trait function and ecological context can facilitate both the detection and the causal understanding of spatiotemporal variation in natural selection.

Transcriptome analyses reveals the dynamic nature of oil accumulation during seed development of Plukenetia volubilis L.

Sacha inchi (Plukenetia volubilis L.) is a shrub native to Amazon rainforests that's of commercial interest as its seeds contain 35-60% edible oil (dry weight). This oil is one of the healthiest vegetable oils due to its high polyunsaturated fatty acid content and favourable ratio of omega-6 to omega-3 fatty acids. De novo transcriptome assembly and comparative analyses were performed on sacha inchi seeds from five stages of seed development in order to identifying genes associated with oil accumulation and fatty acid production. Of 30,189 unigenes that could be annotated in public databases, 20,446 were differentially expressed unigenes. A total of 14 KEGG pathways related to lipid metabolism were found, and 86 unigenes encoding enzymes involved in α-linolenic acid (ALA) biosynthesis were obtained including five unigenes encoding FATA (Unigene0008403), SAD (Unigene0012943), DHLAT (Unigene0014324), α-CT (Unigene0022151) and KAS II (Unigene0024371) that were significantly up-regulated in the final stage of seed development. A total of 66 unigenes encoding key enzymes involved in the synthesis of triacylglycerols (TAGs) were found, along with seven unigenes encoding PDCT (Unigene0000909), LPCAT (Unigene0007846), Oleosin3 (Unigene0010027), PDAT1 (Unigene0016056), GPDH (Unigene0022660), FAD2 (Unigene0037808) and FAD3 (Unigene0044238); these also proved to be up-regulated in the final stage of seed development.

Transcriptome analysis of the medicinally significant plant Fontainea picrosperma (Euphorbiaceae) reveals conserved biosynthetic pathways.

Euphorbiaceae is a large and diverse family of herbs, shrubs and trees that includes a number of species of considerable economic importance as sources of food, medicines and raw materials. One member of this family, Fontainea picrosperma, is the source plant for the diterpene ester tigilanol tiglate, a natural product recently approved as a treatment for canine mast cell tumours. Here we report the development of reference transcriptomes from root and leaf tissues of F. picrosperma, which include core diterpene biosynthesis genes. A total of ~12 Gb of combined clean reads were generated for assembly into 167,566 contigs with a GC (guanine-cytosine) content of ~41%. Gene ontology showed that 2286 and 2504 transcripts were enriched in the cellular process and 2369 and 2529 transcripts were enriched in the metabolic process categories in leaf and root tissue, respectively. The reference transcriptome contains genes coding for core enzymes involved in common secondary metabolite biosynthetic pathways, including the diterpene biosynthesis pathway within the mevalonate (MVA) and 2-C-methyl-D-erythritol 4- phosphate (MEP) pathways. A phylogenetic analysis using these genes found that F. picrosperma clustered most closely to Jatropha curcas. We found a significantly higher concentration of tigilanol tiglate in F. picrosperma root tissue, which correlated with higher levels of gene expression for enzymes associated with the MVA (6 genes) and MEP (7 genes) pathways, and we hypothesise that the initial stages of tigilanol tiglate biosynthesis occur primarily in the roots of F. picrosperma. This study provides a resource for future gene-related biodiscovery investigations in F. picrosperma and diterpene biosynthesis, in particular for tigilanol tiglate and related macrocyclic diterpenes.

Phylogeny and biogeography of Astraea with new insights into the evolutionary history of Crotoneae (Euphorbiaceae).

We investigated species relationships in Astraea, a primarily Neotropical genus of tribe Crotoneae centered in Brazil, using data from the nuclear ribosomal ITS, and the plastid trnL-trnF and psbA-trnH spacers. With all species of Astraea sampled, along with representatives from across Crotoneae, the evolutionary history of Astraea was interpreted in a broader framework, as well as divergence time estimates and reconstructions of ancestral areas and morphological character states for Crotoneae. Our results show that Astraea is monophyletic, consisting of three main clades, and that most of its diversification took place from the Oligocene to the Pliocene, coincident with the formation of the South American "dry diagonal". As for Crotoneae, our data show incongruent phylogenetic positions between the nuclear and chloroplast data for most of its genera, and that the ancestor of the tribe was probably arborescent and might have occupied the Amazon Basin, most likely in moist forest, from which it spread throughout South America in the early Eocene. Ancestral state reconstruction recovered deeply lobed leaves and staminate petals bearing moniliform trichomes as putative synapomorphies for Astraea, whereas the absence or strong reduction of pistillate petals is widespread in Crotoneae and may be a synapomorphy for the tribe.

Larvicide potential of genus CrotonL. (Euphorbiaceae)inbiological control of Aedes aegyptiL. (Diptera: Culicidae)

The reemerging diseases caused by Aedes aegyptiare one of the main public health problems in the world. The control of mosquitoes using larvicidal compounds from products of plant origin is anexcellent alternative. This study aims to evaluate the larvicidal potential of fractions in hexane, chloroform, ethyl acetate and hydromethanol from the ethanolicleaf extract of two species of the genus CrotonL. (Euphorbiaceae) against larval forms of A. aegypti, as an alternative tool to control this vector. Dry leaves of Croton betaceusBaill. and Croton lundianus(Didr.) Müll.Arg. were used for biological tests. The compounds were extracted with ethanol (99.8%). The ethanolic extracts of the leaves were suspended in a methanol / water solution and were successively subjected to the liquid-liquid division process with solvents of different polarities: hexane, chloroform and ethyl acetate, giving rise to the four fractions. Larvicidal tests were performed with the ethanol extract and fractions resulting from the partition. In the study, the crude extract and the fractions showed larvicidal potential, being hexane fractionthe one with greatest activity.Mortality in C. betaceusfractions was up to 40%. Croton lundianuspresented mortality of up to 93.33% of the larvaesubmitted to the test. Data analysis showed larvicidal activity in the crude extract and fractions. The hexane fraction was more effective, especially in C. lundianus.

Genome-wide identification, phylogeny, and expression analysis of the SBP-box gene family in Euphorbiaceae.

BACKGROUND: Euphorbiaceae is one of the largest families of flowering plants. Due to its exceptional growth form diversity and near-cosmopolitan distribution, it has attracted much interest since ancient times. SBP-box (SBP) genes encode plant-specific transcription factors that play critical roles in numerous biological processes, especially flower development. We performed genome-wide identification and characterization of SBP genes from four economically important Euphorbiaceae species. RESULTS: In total, 77 SBP genes were identified in four Euphorbiaceae genomes. The SBP proteins were divided into three length ranges and 10 groups. Group-6 was absent in Arabidopsis thaliana but conserved in Euphorbiaceae. Segmental duplication played the most important role in the expansion processes of Euphorbiaceae SBP genes, and all the duplicated genes were subjected to purify selection. In addition, about two-thirds of the Euphorbiaceae SBP genes are potential targets of miR156, and some miR-regulated SBP genes exhibited high intensity expression and differential expression in different tissues. The expression profiles related to different stress treatments demonstrated broad involvement of Euphorbiaceae SBP genes in response to various abiotic factors and hormonal treatments. CONCLUSIONS: In this study, 77 SBP genes were identified in four Euphorbiaceae species, and their phylogenetic relationships, protein physicochemical characteristics, duplication, tissue and stress response expression, and potential roles in Euphorbiaceae development were studied. This study lays a foundation for further studies of Euphorbiaceae SBP genes, providing valuable information for future functional exploration of Euphorbiaceae SBP genes.

Epigenetic aspects of floral homeotic genes in relation to sexual dimorphism in the dioecious plant Mercurialis annua.

In plants, dioecy characterizes species that carry male and female flowers on separate plants and it occurs in about 6% of angiosperms; however, the molecular mechanisms that underlie dioecy are essentially unknown. The ability for sex-reversal by hormone application raises the hypothesis that the genes required for the expression of both sexes are potentially functional but are regulated by epigenetic means. In this study, proteomic analysis of nuclear proteins isolated from flower buds of females, males, and feminized males of the dioecious plant Mercurialis annua revealed differential expression of nuclear proteins that are implicated in chromatin structure and function, including floral homeotic proteins. Focusing on floral genes, we found that class B genes were mainly expressed in male flowers, while class D genes, as well as SUPERMAN-like genes, were mainly expressed in female flowers. Cytokinin-induced feminization of male plants was associated with down-regulation of male-specific genes concomitantly with up-regulation of female-specific genes. No correlation was found between the expression of class B and D genes and the changes in DNA methylation or chromatin conformation of these genes. Thus, we could not confirm DNA methylation or chromatin conformation of floral genes to be the major determinant regulating sexual dimorphisms. Instead, determination of sex in M. annua might be controlled upstream of floral genes by one or more sex-specific factors that affect hormonal homeostasis. A comprehensive model is proposed for sex-determination in M. annua.

Integrative analysis reveals evolutionary patterns and potential functions of SWEET transporters in Euphorbiaceae.

As a macromolecular substance, sucrose contributes to the plant growth and development. SWEET genes, a group of sugar transporters, are a recently found plant gene family and play important roles in sugar efflux, pollen nutrition, nectar secretion, phloem transport, and seed development. The SWEET genes have been identified and characterized in some plants, but the systematic study in tung tree (Vernicia fordii) was limited. Here, we identified 121 SWEETs in five Euphorbiaceae, and could be divided into four classes with 20 different motifs. Multiple sequence alignment revealed seven transmembrane helixes (TMHs) in the SWEET proteins which were created by an internal duplication of an ancestral three-TMHs unit, connected by TMH4. This study provides direct evidence for the first time for internal duplication in Euphorbiaceae. The large-scale duplication events represented the main driving force for SWEET family expansion in Euphorbiaceae. In addition, we determined the key VfSWEETs for sucrose transport from source to sink tissues in V. fordii and proposed a possible sucrose transport model, which would be helpful for understanding the mechanism of sucrose transport in V. fordii. This study provided a new insight into the evolution, expression and structural variations of SWEETs in V. fordii and four other Euphorbiaceae.

YY males of the dioecious plant Mercurialis annua are fully viable but produce largely infertile pollen.

The suppression of recombination during sex-chromosome evolution is thought to be favoured by linkage between the sex-determining locus and sexually antagonistic loci, and leads to the degeneration of the chromosome restricted to the heterogametic sex. Despite substantial evidence for genetic degeneration at the sequence level, the phenotypic effects of the earliest stages of sex-chromosome evolution are poorly known. Here, we compare the morphology, viability and fertility between XY and YY individuals produced by crossing seed-producing males in the dioecious plant Mercurialis annua, which has young sex chromosomes with limited X-Y sequence divergence. We found no significant difference in viability or vegetative morphology between XY and YY males. However, electron microscopy revealed clear differences in pollen anatomy, and YY males were significantly poorer sires in competition with their XY counterparts. Our study suggests either that the X chromosome is required for full male fertility in M. annua, or that male fertility is sensitive to the dosage of relevant Y-linked genes. We discuss the possibility that the maintenance of male-fertility genes on the X chromosome might have been favoured in recent population expansions that selected for the ability of females to produce pollen in the absence of males.

Early Sex-Chromosome Evolution in the Diploid Dioecious Plant Mercurialis annua.

Suppressed recombination allows divergence between homologous sex chromosomes and the functionality of their genes. Here, we reveal patterns of the earliest stages of sex-chromosome evolution in the diploid dioecious herb Mercurialis annua on the basis of cytological analysis, de novo genome assembly and annotation, genetic mapping, exome resequencing of natural populations, and transcriptome analysis. The genome assembly contained 34,105 expressed genes, of which 10,076 were assigned to linkage groups. Genetic mapping and exome resequencing of individuals across the species range both identified the largest linkage group, LG1, as the sex chromosome. Although the sex chromosomes of M. annua are karyotypically homomorphic, we estimate that about one-third of the Y chromosome, containing 568 transcripts and spanning 22.3 cM in the corresponding female map, has ceased recombining. Nevertheless, we found limited evidence for Y-chromosome degeneration in terms of gene loss and pseudogenization, and most X- and Y-linked genes appear to have diverged in the period subsequent to speciation between M. annua and its sister species M. huetii, which shares the same sex-determining region. Taken together, our results suggest that the M. annua Y chromosome has at least two evolutionary strata: a small old stratum shared with M. huetii, and a more recent larger stratum that is probably unique to M. annua and that stopped recombining ∼1 MYA. Patterns of gene expression within the nonrecombining region are consistent with the idea that sexually antagonistic selection may have played a role in favoring suppressed recombination.